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I - Funciones de KR” en KR”. Límites y continuidad 


El espacio euclídeo R” 


En el conjunto R” = lx, a X, )x, ER k= 1,2,.->,2) definimos la distancia euclídea entre 


dos puntos, x= 0 TN =, y= AAA Y ), así: 


El conjunto R” con esta definición de distancia es el espacio euclídeo R” (en los casos n=1, 
n=2 y n=3 esta definición de distancia coincide con la noción intuitiva). 


La norma de un vector es 


con lo que la distancia de x a y equivale a la norma de su diferencia, d (x, y) = [hc - "| j 


Llamamos bola abierta de centro x€ R” y radio r >0 al conjunto 
Blx,r)= blye R".d(x,y)<rj 


(para n=1 la bola es el intervalo (x A, E + r), para n=2 la bola es un círculo sin incluir su 


frontera o circunferencia de radio r que lo limita, y para n=3 es una esfera, también sin su 
frontera la superficie esférica de radio 7). 


Análogamente, la bola cerrada de centro x€ R” y radio r >0 es el conjunto 


B(x,r)=(y|ye R”, d(x,y)<r) 


(para n=1 será el intervalo [x—r, x+rl, para n=2 y n=3 el círculo y la esfera con su 


frontera, la circunferencia y la superficie esférica respectivamente). 


Decimos que un punto x es interior a un conjunto 4 ER”, o lo que es lo mismo, que 4 es un 
entorno de x, si y sólo si se verifica que existe un real r > 0 tal que B(x,r) C A (es decir, es una 


definición idéntica a la del espacio XX sustituyendo intervalo por bola). Los ejemplos más sencillos 
de entornos de un punto x son las propias bolas B(x,r) . Al conjunto de los puntos interiores de Á 


o 
lo llamamos interior de A, con la notación A. Cuando un conjunto Á coincide con su interior 
decimos que es abierto. Por ejemplo, las bolas abiertas son conjuntos abiertos, las bolas cerradas 
no: los puntos de su frontera no son interiores. 


Dado un conjunto A < /R”, decimos que un punto x de /R”, perteneciente o no a Á, es un punto de 
acumulación de A si y sólo si toda bola abierta B (x,r) contiene algún punto distinto de x que 


pertenece a 4 (o bien, intuitivamente, si existe una sucesión no constante de puntos de Á cuyo límite 
es x). Al conjunto unión de A y el conjunto de sus puntos de acumulación lo llamaremos adherencia 


de A, A. Por ejemplo, la adherencia de la bola abierta B(x,r) es la bola cerrada B (xr) (y el 


interior de ésta es aquélla). Si en R? un conjunto D está limitado por una curva cerrada (o en R? 


por una superficie cerrada), incluyendo o no todos o algunos de los puntos de ésta, el conjunto D 
será la unión de D y los puntos de la curva (o superficie) frontera.. 


Funciones de RR” en R” 
Sea funa función de .R” en RR”, es decir, una función 
f:DCR">R”, 
PA ,%,,.0,%,) = CN 


Para m=1l tenemos una función real, y para m>1 una función vectorial, en ambos casos de n 
variables reales. Si fes una función vectorial, a las m funciones reales de n variables 


Ff/:D>R 
Eo 2,) =P, 


para ¡= 1,2,...,m, se las llama funciones coordenadas de f. 


En las funciones reales, para n>2 el aumento del valor de n ya no introduce diferencias 
importantes de concepto; pero en el caso particular n = 2 la gráfica de /; es decir, el conjunto 


(E y,z)((, y)eDc R?z= (x, ») está en R? y es intuible, dibujable. Pueden hacerse 


representaciones mediante cualquier sistema, por ejemplo utilizando secciones por planos: en 
particular los coordenados, o bien planos z =k, secciones éstas últimas que proyectadas sobre el 
plano z = 0 constituyen las llamadas curvas de nivel. 


Ejemplos. Las funciones definidas por las expresiones 


xy 


faoy=x+y, floy)=xy, f(x, y) =maxdlxa, lvl), £(0,y)= pe 


cuyas gráficas son como indican las figuras, 


0 
a 


TN 
SALA ART 
RAR 7 
AARRATA yy 


SS 


Idea de curvas y superficies: ecuaciones paramétricas 


Entre las funciones vectoriales de varias variables hay tres casos especialmente importantes desde el 
punto de vista geométrico. 


A pe 3 
Siendo 7 un intervalo de la recta real, supongamos una función f:I=>R', 


Fs (4,0, AO) = (x(0), y(0)). La imagen (¡no la gráfica!) de f, 


((x(0)y(0)|r e DER”, 


es, si ftiene una cierta regularidad (como pasaba en una variable, y según irá quedando claro más 
adelante), una curva plana. A las expresiones 


10 >. e x(1) 
o simplemente 
pq 


se las llama ecuaciones paramétricas de la curva. 
Ejemplos 
La función 
f:[0,277) >R” 
lb (cost, sen 1) 


tiene como imagen la circunferencia de centro el origen y radio 1; o con otro lenguaje, las 
ecuaciones paramétricas de ésta son: 


fx =COSÍ 


| y =sent 


16 [0,277). 


Asimismo, las ecuaciones 


x=r(t- sent) 
leR, 
y =r(1-cost) 


son las de la curva llamada cicloide, descrita por un punto de una circunferencia de radio r que 
rueda sobre el eje de abscisas. 


Supongamos ahora, siendo / un intervalo de la recta real, una función 
PASR, FO= (40. (0,4(0)=(x(0),y(0),2(0)); 


su imagen será una curva en el espacio de tres dimensiones, una 
curva alabeada. 


Ejemplo. La hélice: 


x=acosÍ 


y=asent tek. 


A e iaa 3 
Finalmente, si fes una función de /XJ en R”, donde 1 y J son 
intervalos, 


Fu, y) = (x(u, v) ylu, v) zlu, ), 


o bien en paramétricas, 


=x(u,v) 
y=y(uv) (uv)je1IxJ, 
z=z(u,v) 


tendremos una superficie. 
Ejemplos 


¿Qué superficie es la de ecuaciones paramétricas 


x =C0SUSCNV 
y=senusenv (u,v)e[0,27)x[0,71)? 


Z=COSV 


% O 
(calcúlese x?* + y* +20). 
¿Cuáles son las curvas de nivel z=4, z=1l y z=0 de la superficie 


x= (1 +u)cos v 
y=(1-u)senv  (u,v)e RxÍ0,277)? 
z=g" 


Las figuras dan idea de cómo es la superficie (para -2<u <2) y cuáles son las curvas pedidas. 


E 
Se 
7) 
/ 
Hi] 
Y 
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SS 


Límites 


Es la noción básica del cálculo diferencial en varias variables, como lo era del cálculo en una sola. 


Empezaremos por el caso de funciones de R? en R, funciones reales de dos variables reales, 
fácilmente generalizable a funciones reales de n variables. Y consideraremos funciones definidas en 
un dominio D abierto y conexo (un conjunto es conexo si y sólo si dos puntos cualesquiera del 
mismo pueden unirse mediante una curva incluida en él; o más formalmente, si dados dos puntos 


cualesquiera del conjunto D existe una aplicación continua (Y de un intervalo [a,b] CR en D tal 


que dichos puntos sean precisamente pla) y plb)*. 


Definición. Sea f : D=>R, donde D es un dominio de R?, (a,b)e D (adherencia de D), y 
l € R.Se dice que les el límite de fen (a,b) , con notación 


ea Fasj=l, 


(x,y)> (a,b) 
si y sólo si para todo intervalo (bola abierta) U= (1 —E,l +€) existe una bola abierta 
V= B((a,b),9) tal que, si un punto (x, y) distinto de (a,b) pertenece a Y, entonces F(x,y) 
pertenece a U, O lo que es lo mismo: 


lin 184 )]=] 
(o e (x, y) 


si y sólo si, para todo £>0 existe un 0 >0 tal que, si 0<|(x, y) -(a,b)||< 0 entonces 


(uy) -l|<e. 


En lo sucesivo, y para aligerar la notación, en lugar de Le Ay) escribiremos simplemente 
(uy (a,b 


limf (x, y) , expresión que no crea ninguna ambigúedad. 
(a,h)' ¿ 
La generalización de esta definición para n>2 es inmediata. Y para el caso general de funciones 


vectoriales (es decir, en R”) y con la nueva notación tendremos: 


Definición. Sea f:D=>R"” donde D es un dominio de R”", a=(a,a,,-,a,Je D y 
sl 1 


Je R”.. Se dice que / es el límite de fen a, con notación 


m 


limf(x)=1, obien limf(x%,...,x,)=(l,,..l,,), 


a 


si y sólo si para toda bola abierta U = B(1,€) en R” existe una bola abierta Y = B(a,9) en R” 


* En rigor ésta es la noción de conexo por arcos, más restrictiva que la de conexo. Pero, teniendo en cuenta que en 
abiertos ambas nociones coinciden, es suficiente a nuestro nivel. 


tal que, si un punto x distinto de a pertenece a V, entonces F(x) pertenece a U, O análogamente: 


limf(x)=1, o bien lim f(x,,...,x,)=(l,...,L,,), 


Fx) <e. 


si y sólo si para todo € > 0 existe Ó > 0 tal que O< llx—al| <óÓ implica | 


Si 17 (x) — Il|< E, como para ¡=1,2,...,m es 


m 


401) =|1()-1<0, 


resulta 
EA AAA O lim f, ael, para 1,3, 


la propiedad recíproca, es decir, 


se demuestra también sin mayor dificultad, con lo que se obtiene: 


el límite de fen un punto a es / si y sólo si el límite de cada función coordenada f, es la 


correspondiente coordenada ? del límite. 


Operaciones con límites 


Todo lo visto en funciones reales de una variable en relación con suma, producto, valor absoluto, 
función inversa (para el producto) y exponenciación es aplicable a las funciones reales de varias 


0 

, 2 ñ y A 0 : s a 
variables, resultando también los casos de indeterminación 0 y 0. Y si, de la forma obvia, 
ampliamos como allí el concepto de límite a límites infinitos, aparecerán también los casos co—0o, 


E 0 co 
o0:0,—,0 y I”, 


oo 


Límites según una trayectoria. Límites radiales 


Para el cálculo de límites en los casos de indeterminación suele ser útil el concepto de límite según 
un conjunto: su definición es análoga a la de límite, pero considerando sólo, en los entornos del 
punto a, los puntos x que pertenecen a un cierto conjunto E, del cual a es punto de acumulación. En 
particular consideraremos los casos en que E es una curva (límites según una trayectoria), 


posiblemente una recta (límites radiales). Evidentemente, si lim f(x) = 1, entonces / es también el 


x=ad 


límite de fen a según cualquier conjunto. En particular, la existencia e igualdad de todos los límites 


radiales es condición necesaria (pero no suficiente, como pronto veremos) para la existencia de 
límite, y puede ser útil a veces para concluir rápidamente que una función no tiene límite. 


Ejemplo 1. Sea 


¿tendrá límite en el punto (0,0) 2? Para cada dirección y = mx resulta 


y 
, . mx? m 

lim f (x,mx) =lim—_—== == » 

x>0 x>0 (1 +m de l+m? 


es decir, todos los límites radiales existen pero no son iguales, varían con m. Luego la función no 
tiene límite en ese punto (véase gráfica y curvas de nivel en la página 3). 


Ejemplo 2. Sea ahora 


2x y 
f(x,y)= z 
ur) = 
24x a ; ; 
Puesto que f(x, Ax) ==, los límites radiales en (0,0) existen y son todos iguales a 0; pero 
q ¿ y? e P 


: : : 2 
si consideramos las trayectorias y = 4 x” resulta: 


4 
lim f (xx?) =lim_ 422 
x>0 x>0 x +u Y 14 ul" 


los límites existen para todas estas trayectorias parabólicas, pero no son 
Ha iguales, dependen de 4: la función no tiene límite en (0,0). Es decir, la 
condición necesaria recién vista no es condición suficiente. 


Cambio de coordenadas 


El límite de una función en un punto puede expresarse en coordenadas polares combinadas con una 
traslación al punto en cuestión. Recordemos que en estas coordenadas el punto (x, y) viene dado 


por el par (r, 0) , con las relaciones 


x=rcos0 


y=rsen0. 


La bola abierta B((0,0), 5) se convierte en el conjunto [(7,0)| r00S0< 27), con lo que 


e Fx, y) se convierte en lim f (rcos0, rsenó); la existencia de éste con independencia de la 
0,0 r>0 


libre variación de O equivaldrá a la existencia de q Fl y) con el mismo valor. 
0,0 p 


Ejemplo 1 
$. 3 $ ad 
E O a ES A 
lim 52 = A = limr (cos' 0O+sen 0) =0, 
(0,0) x?* + y” r>0 r r>0 
ya que la función del paréntesis está acotada. Si el límite fuera en un punto (a,b) habría que 


combinar el paso a polares con una traslación a este punto, es decir, hacer el cambio 


x=a+rcos0 


y=b+rsen0 
y el límite seguiría siendo para r >0. 


Ejemplo 2 


/ x+ o =Y+rcosO+YV+rsend0 . cosO+senO 

lim ar A de E O mM —_—_—_—_—_—_———_ 

Ea -y+1 r>0-Y+rcos0-Y-—rsen0+1l ">0cos0—senO 
y 


que evidentemente no existe. 


Si se tratara de una función real de tres o más variables, todo lo dicho sobre límites según un 
conjunto seguiría siendo válido. Pero para n= 3 puede usarse también un cambio de coordenadas. 
Hay dos sistemas muy usuales, y que en todo caso tendrán aplicación más adelante, que son las 
coordenadas cilíndricas y las coordenadas esféricas. 


En las primeras, un punto (x, y,z) se representa como (r,0,2), 
siendo 

x=rcos0 

y=rsen0 
y quedando la z sin cambios. Es decir, consisten simplemente en 


representar en polares la proyección sobre el plano OXY, (x, y). (La 
ecuación r = cte. es la de un cilindro de revolución de eje OZ.) 


En las coordenadas esféricas el punto (x, y,z) se representa como ( P,0,9), donde p es 
ly, Dl, 0 el ángulo formado por el vector En y,0) con el eje OX y Q el formado por el vector 


(x, y, z) con la parte positiva del eje OZ, medido a partir de ésta. Es decir, 


x= pcosOsenq 
y =psenOseng 
Z=pcosq 


(la ecuación P =cte. es ahora la de una esfera de centro el origen.) 
Ejemplo. La función 


e al 
(x, y,z) ER 2 2 y > 
CEA 
¿para qué valores de los reales positivos p, q y s tiene límite en (0,0,0)? Pasando a coordenadas 
esféricas y con un razonamiento semejante al caso de las polares en el plano será, dado que 
2 2 2 2 
o O O 
Ese Ñ 
. P s-2 
lim 3 = lim p*"""9(0,p), 
(0.00) + y 42? p>0 


donde D(0, p) es un producto de senos y cosenos de O y de (, es decir, una función sin límite 


pero acotada, luego el límite buscado existe y es O si y sólo sí p+q+s>2. 


Continuidad 


El concepto de continuidad es aquí análogo al que teníamos en el caso de las funciones reales de 
una variable. 


Definición. Una función es continua en un punto si y sólo si su límite en él existe y es igual a su 
valor. Y es continua en un dominio D si y sólo si lo es en cada uno de los puntos de éste. 


Recordando las propiedades de los límites resulta inmediatamente: 


una función de R” en R” es continua si y sólo si lo son todas sus funciones coordenadas; la suma 
y el producto de dos funciones reales continuas, así como el valor absoluto, la inversa (para el 
producto) y la exponencial de una función real continua son funciones continuas. 


Ejemplos 
Acad Xy ; / 
La función (530) += > €s continua en todo el plano real excepto en el origen de 
Ry 
coordenadas, en el que no tiene límite. 
il FEprl De ) , 
La función (3, y) + — 7 €s también continua en todo el plano real excepto en el origen de 
Xx +y 


coordenadas, en el que no tiene límite finito. 
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E 
6 xyz , ñ , , 
La función ( o y,z) > E que es la vista más arriba en el caso particular 

e 
p=q=s=2, es continua en R*—[0,0,0), pero si completamos su definición haciendo 
F(0, 0,0) =0 (prolongación por continuidad) tendremos una función continua en todo el espacio 
R? 


Composición de funciones continuas 


Teorema. Sean f:D>R” y g:E=>R”, con D y E dominios respectivos de R” y R”, 
F(D)< E, f continua en D y g continua en £. Entonces la función compuesta go f:D=>R” 


es continua en D. 


DecP Ec Rr”" sr 


Dor” SISI > A 


Demostración. Análoga a la que teníamos para funciones reales de una variable. Sea un punto 
genérico Y =(x,,x,,""",X,)€ D, y sea F(X)=Y=(y,),""",),)€ E; para toda bola abierta 
W = Blg (y). e) == Blg o f(x),€) CR” existe, por ser g continua, una bola abierta 
V= B(5, 7) CR” tal que ye V => gl») EW ; pero para esta bola Y existe, por ser f continua, 
una bola U=B(X,9) R” tal que xeU => f(x)eV ; luego en definitiva, para toda bola 
abierta W= Blgo f(x). e) CR” existe una bola abierta U=B(x,9)=R tal que 


xeU= go F(x) EW,y gof es continua en x,es decir, en D. M 


Imagen de una función continua 


De forma semejante a lo que ocurría en funciones reales de una variable, se demuestran también las 
dos propiedades que siguen. 


Teorema. Si una función es continua, la imagen de un conjunto K cerrado y acotado (también 
llamado compacto) es un conjunto cerrado y acotado (llamamos acotado a un conjunto de R” tal 
que las distancias de todos sus puntos al origen de coordenadas forman un conjunto de números 
reales superiormente acotado). 


Teorema. Si una función es continua, la imagen de un conjunto conexo es un conjunto conexo; si se 
trata de una función real, por lo tanto, la imagen será un intervalo (o, caso trivial, un conjunto de un 
solo punto si la función es constante). 
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Il — Derivadas y diferenciales 


Derivadas parciales de una función real en un punto 


Sea funa función real de dos variables definida en un dominio D, y La Yo) un punto de D. Si 
fijamos el valor de una de las variables, por ejemplo y = y, , queda definida una función de una 
variable, p(x) = fl, Y» pS en un entorno del punto x= x,; y análogamente, haciendo x= x, 
queda definida w(y)= FUE en un entorno de y= y,. Pues bien, si estas funciones son 


derivables en esos puntos, sus derivadas son las derivadas parciales de fen (x, Da ). 


Definición. Sea f'una función real definida en un entorno de un punto (xo, Hd, Decimos que f 
admite derivada parcial respecto de x (o que es derivable respecto de x) en (Ms, ) si y Sólo si 


la función de » 


Fla +h, Yo ) Dé Fs Yo ) 
h 


definida para h en un entorno de 0, tiene límite finito para h>0; a dicho límite lo llamamos 
entonces derivada parcial de fen pd le con la notación f, (tod ] o bien > or Jo): 
PE +h, Yo) f(%o, Yo) 


F¿(%o Yo) ss lim h ; 


De forma totalmente análoga decimos que fadmite derivada parcial respecto de y en (x, R Yo) si 


y sólo si existe 


k =— , 
1x9) =pin ot Poyo. 


número real que es la derivada parcial de frespecto de y en Ea ¿Yo ¡3 


Estas derivadas miden, naturalmente, la variación de fa lo largo de las direcciones Y = Py y 
X=X,. Y su existencia implica la continuidad de las aplicaciones P y YW de las que 


hablábamos, pero ésta no implica la continuidad de f. 


Ejemplo 1. Sea otra vez 
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Estudiemos sus posibles derivadas parciales en el origen: 


00-92 20-200) =0, e007= 9 OH LO _ pp 


h>0 h 10 k>0 k k>0 k 


admite derivadas parciales en el origen, aunque como vimos no es continua en él. 


Ejemplo 2. Sea ahora 


1(0,1)= 422 si (xy) (0.0) 
x+y” 
f(0,0)=0 
Dj F(n,0)- (0,0). h e LOE)= FOO] 1 
A a AA a a 


que no existe; fadmite en (0,0) derivada parcial respecto de x pero no respecto de y. 


Todo lo anterior se generaliza sin dificultad a funciones de R” en R. Podrán existir ahora hasta 
n derivadas parciales, una para cada una de las » variables. 


Ejemplo 3. Sea la función de R*enR 


f(x, y,z)= yx? + y? +z?. 


En el origen, 


h,0,0)- f 
£.(0,0,0)=tim2 220.0) 4(0,0.0) _ y, Al. 
ñ h>0 h h>0 Y 


que no existe; f no es derivable respecto de x en el origen, y tampoco respecto de y ni de z. 
Interpretación geométrica en el caso de dos variables 


De la definición se desprende que la derivada 
F.(x, yo) de una función real de dos variables es 
la pendiente de la recta tangente en el punto 
lo Yo) a la intersección (curva) de la gráfica de f 
(superficie) con el plano y = y, (ver figura). Y lo 


mismo con la otra derivada Fit y el plano 


Xx =Xp. 
Derivadas direccionales 


Lo mismo que hemos definido, a partir de f; dos funciones de una sola variable fijando el valor 
de x o de y, podemos hacerlo ligando las variables x e y de forma que el punto (x, y) recorra 
una semirrecta engendrada por el punto (xo, Yo) y el vector unitario v= (cos 0,sen0) , CUYOS 


puntos serán Ls +rcos 0, y, +rsen 0), con r positivo; queda así definida, para r >0, una 
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función de una sola variable r, x(r)= f(x, +rc0s0, y, +rsen0), cuya derivada en r=0, 


x'(0)= tim 2LI=210) 


, es, si existe, la derivada de fen (Xo, Yo) en la dirección y sentido 
r>0+ Y 


del vector y. 


Definición. Sea f una función definida en un entorno de un punto (x, Yo ). Llamamos derivada de f 


en la dirección 9 en (x,,,), f, (%.,y,) (6 f, (%o y). siendo v = (cosO,sen0)), a 


lim FS, +rcos0, y, +rsenO)- f(x, yo) 


r>0 r 


si dicho límite existe y es finito. 


Esta derivada mide la variación de fcuando el punto 
(x, y) varía desde (to, Yo) en línea recta y en el 
sentido O, y su existencia equivale a la existencia de 
recta tangente a la gráfica de f' según la dirección O, 
recta cuya pendiente será el valor de f, (e Yo) N 


La definición se generaliza sin dificultad para más 
variables. La dirección vendrá dada ahora por un 
vector unitario v=(v,v,,""",,), y la derivada 


direccional en un punto a = la,,a,,---,a,) será 


r>0 vd 
y su significado será el mismo: medirá la variación de fen a en la dirección de v. 
Función derivada parcial 


Lo mismo que en el caso de funciones de una variable, si una función fadmite derivada parcial, 
por ejemplo respecto de x, en todos los puntos de un dominio D, queda definida en D la función 
derivada parcial, 


F:D>OR 
(y) > £. (xy), 


y lo mismo para la variable y. 


En los ejemplos de más arriba hemos estudiado las derivadas en puntos en los que, por la 
particular definición de las funciones, teníamos que aplicar directamente la definición. Pero en 
la mayoría de los casos ello no es necesario, basta aplicar las reglas habituales de derivación; 
recordemos que una derivada parcial no es más que la derivada de la función de una variable 
que se obtiene de la dada suponiendo constantes todas las variables excepto una. 


Así, en los ejemplos 1, 2 y 3 de las páginas 12 y 13 tendremos respectivamente, para puntos 
distintos del origen, 


E 2 3 2 2 3 0 2 
EA A o E, 
(1? + y?) (x? + y?) (1? + y?) 
3 + y) 2 (x + y?) x(x*+3xy?-2y?) 
(Peg (+y2) 
2er) 2 (+) 2 (13), 
a en | 


qm +y +2 Ja +y +2 :b57,2)= et +y+2 


Caso de funciones vectoriales: matriz jacobiana 


f,(x,y)= 


Lap 2)= 


Si se trata de una función vectorial, por ejemplo (tomando n= 2 por simplificar notaciones) 


F : R? Y R"” 
Fc y)=(4 (oy), £ boy): f, by), 


las definiciones de derivabilidad y de derivada parcial serán las mismas. Tendremos ahora, en 
un punto genérico ¡YA y), 
: e ARA] IS 
Ú h=>0 h 
7 


y A Ln) elas) 
1>0 h 


iS E ESIA SIE : des (x, y)) > 
y análogamente 


fia y)= (Ly (3), Loy (0,9) Lay (2, 9)) 


Las derivadas parciales en un punto son vectores de R”, y las funciones derivada parcial, 
> ; 2 $ % ¿ 
funciones de R* en R” como f. Aparece entonces el concepto de matriz de derivadas, o 
matriz jacobiana, 


Df (x,y)= 
pe (x, y) Por (x, y) 
matriz mX2 en este caso, mXn en el caso general de una función de R” en R”, concepto 


que desempeña un papel muy importante en el cálculo diferencial (e integral) en varias 
variables, como iremos viendo. 
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Derivadas parciales de orden superior 


Lo mismo que en una variable, las funciones derivada parcial pueden a su vez ser derivables, 
con lo que nos encontramos con el concepto de derivada parcial segunda. En dos variables 


tenemos: 


94 


fi. 37 ): derivada parcial segunda respecto de x dos veces; 
2 
): derivada parcial segunda respecto de x y de y; 


0 
F,, (6 
*  0xdy 
2 
Pos (ó 3 — ): derivada parcial segunda respecto de y y de x; 
Jydx 


o) 


¿e (6 5 ): derivada parcial segunda respecto de y dos veces. 


dy? 


Ejemplo. Sea de nuevo 


Fey) => si (x, y) + (0,0) 
x"+y 


/(0.0)=0 


para la que ya hemos calculado, fuera del origen (0,0), 


' xy", Xx” —xy 
¿ano aa 
(ax + y ) (x + y?) 


derivando de nuevo, 


dal (a+ y) (1? + y) | 
Xx (1 + y?) (+ y?) 


(3-9 e+ y] a+ yema) tq 6xy? y 
AAA A 
(x + y ) (x + y ) 


2xy (+ y) ay (+ 9?) (0 y?) _ 6 y 4 2ay' 


f, (x,y)= SSA " 3 
(1? + y?) (3? +3?) 


Permutabilidad de las derivaciones 


Vemos en este ejemplo que las derivadas f,, y f,, resultan ser idénticas; ello no es 


casualidad, como expresa el siguiente 
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Teorema. Sea f una función definida en un entorno U de un punto (x,, Yo) y tal que existen en 
Uf. f, y f,, [respectivamente f,,], siendo ésta última continua en om Pal. Entonces 


[respectivamente f 


existe en ese punto la otra derivada mixta f, w ] y se verifica que 


MEN SS Ea (0: Yo) . 


Omitiremos la demostración. 


Es decir, según este teorema, la posibilidad de que f,, 


y Fs supuesto que ambas existan en 


un punto, tengan en él valores distintos exige que ambas sean discontinuas en dicho punto. 


Desde luego este resultado se generaliza al caso de funciones de » variables, puesto que si, por 


ejemplo, tenemos tres variables, las derivadas f,, (%o, Yp,Zp) y Felt Mis) no son sino las 


derivadas mixtas de una misma función de dos variables, Fx, EA 7 


De la misma forma que las derivadas parciales segundas se definen las terceras, cuartas, etc. Y 
a ellas se aplica también la propiedad de ser indiferente el orden de derivación (supuesto que se 


cumpla la condición de continuidad). Por ejemplo, las derivadas f,,, (rd) Y La EM Yo) no 
son sino las derivadas (4), (oo) y (f. ), (o Yo) , es decir, las derivadas mixtas de una 
misma función, f,. Y lo mismo para más de dos variables. Si Fes, por ejemplo, una función de 
tres variables, las derivadas terceras La (xo, My Za) y Ll tl como derivadas 
(£.),. (o: Yo:Z0) y (f, ),, (%: Yo» Zo) de una misma función de dos variables, f,(x,,y,2), 
están dentro del alcance del teorema y son iguales. 


Si la condición de continuidad se satisface, pues, una función de tres variables tendrá sólo seis 
derivadas parciales segundas, las de subíndices xx, yy, 2Z, XY, YZ y Zx, diez derivadas terceras 


(¿cuáles?), etc. 
Diferenciabilidad de una función real de varias variables 


Llegamos al concepto básico del cálculo diferencial, que es la verdadera generalización del 
concepto de derivabilidad en una variable (recordemos que allí derivabilidad y 
diferenciabilidad eran conceptos equivalentes). 


Definición. Sea f una función real de dos variables definida en un entorno U de un punto 
(x, ; dl: Decimos que fes diferenciable en (Lo. Fo) si y sólo sí existe una aplicación lineal 


A:RESR 
A(h,k)= ah + bk 


tal que se verifica que 
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(obsérvese que se trata del límite de una función de las variables h y k exclusivamente, es decir 
del límite para (h,k) => (0,0) ). 


Llamando 


F(%o +Hh, Yo HK)- (xo, yo) -A(M,k) (h k) 


la diferenciabilidad puede expresarse también diciendo que existen dos números reales a y hb 
tales que se verifica la igualdad 


> Fl +hy +k)- f(x, 1) =ah+bk+/n? +? e(h,ko), 
donde € es una función que verifica 


lime(h,k)=0. 

(0.0) 
La definición sería análoga para un número mayor de variables. Para tres, por ejemplo, la 
aplicación lineal A se definiría A(M,k,1)=ah+bk+cl, etc. 
Está claro que, como en el caso de una variable, la diferenciabilidad de una función en un punto 
expresa la existencia para ella de una buena aproximación lineal local, en el sentido de que la 
diferencia de los valores de f y de la aplicación lineal en puntos de un entorno del dado es un 
infinitésimo de orden superior a la distancia entre aquéllos y éste. (Aproximación que 
geométricamente es la ecuación del plano tangente a la gráfica de f en E ¡Em Y)) , de 


forma análoga a lo que ocurría con la recta tangente en las funciones de una variable; 
volveremos sobre esto.) 


Propiedades: continuidad y derivabilidad 
Teorema. Si una función f es diferenciable en un punto, entonces es continua en él. 


Demostración. Tomando límites en la expresión (1) resulta 
lim PQ +A yo +Kk)- flo, yo)= limfa +54 +4 h? + e(h,k))= 0.m 


Teorema. Si una función f es diferenciable en un punto, entonces admite en él derivadas 
parciales. 


Demostración. Haciendo k = 0 tenemos 
Fo +h, yo)= (xo: Yo)= ah +|ne(1,0), 
y operando y tomando límites, 


£. (xo y.) =timL Lotto 90)= 4019) - 


h>30 h h>50 
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es decir, existe la derivada parcial respecto de x y su valor es precisamente el coeficiente a. Y 


lo mismo se comprueba para la otra derivada parcial: f, (Mosa) =b.M 


Análogamente para mayor número de variables. 
Diferencial de una función real 


A la aplicación lineal a la que hemos venido llamando 4 se la llama por definición diferencial 
de fen el punto considerado. En el caso de dos variables, llamando dx y dy a lo que hemos 
llamado hasta ahora h y k, y dfa A, tenemos 


df (Sip Yo ) =hf, Em Y ) $ K£, En Yo ) 
como expresión habitual de la diferencial, o incluso, si llamamos f(x, y) ==, 


dz EN Yo) a hf, [ys Yo ) 3 Kf, e Yo) . 
Para tres variables tendríamos 
df (Xo, Yo, Ea) =hf, e Yo % ) e (ds Zo ) +/f. Zo) ; 
y análogamente para n variables. 
Ejemplo 1. Sea la función de R? en R 


+y 
NM 


F(0,0)=0 


si (x,y) (0,0) 


Estudiemos su diferenciabilidad en el origen. 


17 
ll : : 
h,k)- f(0,0)-Af.(0,0)—kf (0,0 
p h>0 h 3 3 (0,0) |? +? 
TA: 
2,1. 343 3.3... hk*-Hók_, -risenOcos0(sen0+cos0) 
.n h? a = a An 6 3 ES lim 37 ; 
yn 3(1 +1?) 


que no existe: fno es diferenciable en (0,0) (la gráfica de fno tiene plano tangente en el origen 


de coordenadas). 


Ejemplo 2. Sea ahora 
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h* 
—-—(0 ] 
h,k)-f (0,0) -kf. (0,0) -Xxf (0,0 
1,(0,0)=timl=—=0, f,(0,0)=0; y E) 0.0) (0.0) (0,0) _ 
h=>0 (0,0 he + fe? 
h*+k* 
2 2 4 4 Aloryt 4 
tim BELLE — iy PARE — pi sen 8408 0) _ y, 


LN PD O je r>0 pi 


esta función sí es diferenciable en el origen, y su diferencial en él es dz =0 (o bien, la gráfica de f 


tiene plano tangente en (0, 0,0) : el plano z=0). 


Funciones de clase C” 


Lo mismo que en el caso de una variable, se llama así a las funciones que tienen derivadas parciales 
continuas hasta las de orden » inclusive. 


Hemos visto que toda función diferenciable es continua y admite derivadas parciales de primer 
orden. Sin embargo, éstas no son necesariamente continuas. Sin embargo sí es cierto el recíproco. 


" . n ., 1 E A) E 
Teorema. Si una función es de clase C” en un punto, entonces es diferenciable en él. 


Es el modo más sencillo de comprobar la diferenciabilidad de una función en muchos casos. Así, las 
3 4 4 
sb a d : e ; xo +y x +y 
funciones de los ejemplos 1 y 2 de más arriba, definidas por las expresiones === Y ==; 
A A e 
respectivamente, tienen, fuera del origen, derivadas racionales cuyo denominador no se anula, es 
decir, continuas, luego ambas funciones son diferenciables en todos esos puntos. 


$ A 3 ., . . I 2 
Como ejemplo, en cambio, de función diferenciable en un punto pero no de clase C” en él tenemos 


F(xy)= (x+ y) sen 


la función X+ Y en cualquier punto de la forma (Xy, Xy ) , es decir, de la 
Ff(x,-x)=0 
recta x+ y =0: compruébese, particularizando por ejemplo en el punto (1-1), primero, que la 


función es diferenciable en él, segundo, que las derivadas parciales existen en un entorno del 
mismo, y tercero, que no son continuas (no tienen límite) en él. 


Derivadas direccionales de una función diferenciable 


Sea funa función diferenciable en un punto (x, ¿Yo ). Estudiemos su derivada en ese punto según la 
dirección y sentido dados por el ángulo O con la parte positiva del eje OX, es decir, por el vector 
unitario v=(cos O, senO). Será, si existe, 


A f(x, +rcos0, y, +rsenO)- f (xy, Yo) 
JO 0-40) — i 


r>30 r 
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Pero, por la diferenciabilidad de /, haciendo h=rcos0, k=rsen0, con lo que Vh*+k* =r, 
elh, k) = elr,0) con lim e(r,0) = 0, podemos escribir 


ES +rcos0, y, +rsen0)- F(%0 Y) =f. (xo, Yo)r cos0 + f, (x, y, )r senO +r e(r,0), 
y dividiendo por r y tomando límites, 


fo (o, Yo) Ss F. (xo Yo) cos 0 + f,(Xo, Yo) sen0 WA ENANA v, 


el producto escalar del vector (£.(%:Y0)..£, (%o, Y.) por el vector unitario v; una función 


diferenciable tiene derivadas según todas las direcciones y sentidos, dadas, en función de las 
derivadas parciales, por la expresión que acabamos de ver. Y lo mismo para mayor número de 
variables, con la dirección dada por el vector unitario v= ¿A v,), Il =1. Para n=3, por 


ejemplo, 
Á (Xo» Yo» Zo) =f, EEN + Li (Ao, Y, ¿Zo )V» EÑ LAA 


Gradiente 


Definición. Dada una función f de dos variables, diferenciable en un punto, llamamos gradiente de f 
en dicho punto al vector 


CACA 


(léase “Vf” como “nabla de f”). Y lo mismo para n variables. El gradiente es, pues, Df, la matriz 
(1 eS 2) (ó (1 X n) ) de derivadas, considerada como vector del espacio euclídeo RR" de 
Según lo que acabamos de ver, la derivada direccional es el producto escalar del gradiente por el 


vector unitario en la dirección y sentido considerados (se la llama también derivada según un vector 
unitario): 


fol o Yo) =VF (o Yo) (cosO,senO) =Vf (Xp, Yo): v, 
o bien, para cualquier número de variables y sin especificar el punto, 


JJ y da Vf ió v 6 

Recordando que Vf -v= [V7] : Il -COSA= INZA| -cosQr, donde (% es el ángulo que forman Vf y 
v, tenemos una interpretación geométrica para el gradiente. La derivada f, =Vf -v, que mide la 
variación de f en la dirección de v, es máxima cuando cos =1, es decir, cuando 4=0 y la 
dirección y sentido de Vf y v coinciden; luego el gradiente señala la dirección y sentido del 
máximo crecimiento de f. 


Interpretación geométrica en el caso de dos variables 


Como ya hemos dicho, la diferencial, dz(x,,y,)= f, (xo, y) dx + f, (xo, Y) dy , es en este caso 


la ecuación de un plano referida a unos ejes (dx, dy, dz) paralelos a los dados y que pasan por el 
punto em Yo» Z(Xo, yn ), plano que es, por definición, el plano tangente a la gráfica de fen dicho 


punto (obsérvese que contiene las dos rectas tangentes dz= f, (x,, Yo) dx y de=f,(%, Yo) dy, 
paralelas a los planos OXZ y OYZ): la 
diferenciabilidad equivale a la existencia de 
plano tangente. Pasando a los ejes 
originales de coordenadas, la ecuación de 
dicho plano será 


OD 2-2%=2,(x-x)+z, (yy). 
La existencia de plano tangente se 
traduce en que, primero, existen todas 
las rectas tangentes a la superficie en el 
punto, es decir, todas las derivadas 
direccionales, y segundo, en que 
aquéllas son coplanarias, es decir, las 
derivadas direccionales están ligadas 
(por la expresión hallada) con las 
derivadas parciales. 


En cuanto al gradiente, al señalar la 
dirección y sentido del máximo 
crecimiento de f señala también la 
dirección y sentido de la recta de 
máxima pendiente en el plano tangente, 
y será por ello ortogonal tanto a la recta 
de pendiente nula en dicho plano como a 
la correspondiente curva de nivel, puesto 
que aquélla es tangente a ésta. 


Ejemplo 1. La figura representa la gráfica de la 
función del ejemplo 2 de la página 19, 


S 
dl 
Y 


RA SS Je P 
DOSIS 4 f(x, y)=5 2 si (x, y) + (0,0) 
E Y 
1(0,0)=0, 


junto con el plano tangente en el origen, z=0; es obvio que todas las rectas del plano que pasan 
por el origen son tangentes a la superficie. 


Ejemplo 2. Aquí tenemos en cambio la función del ejemplo 1 de la misma página, 


3 
; PE a 
flay)= a - si (x,y)+(0,0) 
x +y 
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y el plano de la forma (1) de la página 22, 


x= 
g= 


, que no es tangente en este 


caso: hay rectas tangentes a la superficie 
que no pertenecen a él. 


Ejemplo 3. Sea la función f de R? en R 
2 rd, definida por la expresión 


Ay)=*-39y + y +7x+6; 


en el punto (5,p)= (1,2) el gradiente de f es 
Vf(1,2)=(=1L1), el valor de la función f(1,2)=10, la 
curva de nivel para este valor de z es la hipérbola 
=x?-3xy+y?+7x-4=0, el plano tangente el de 
ecuación z—10=-(x-1)+ y—2 y su intersección con el 
plano z=10 la recta =(x-1)+y-2=0, x-y+1=0. 
En la figura se ve representado todo esto en el plano 
z= FLZ) =]0, observándose que la recta es tangente a 


la curva y que el gradiente es ortogonal a ambas. 


Diferenciabilidad de una función vectorial de varias variables 


Si se trata de una función vectorial, f : R" => R”, la definición de diferenciabilidad es 
formalmente idéntica a la que hemos visto para funciones reales. Dado un punto 
A x+h=(x,,x,..,x,)+ (M,,h,...h, ) variable en un entorno del 


primero, decimos que f es diferenciable en x si y sólo si existe una aplicación lineal 
A:R” => R” tal que 


pim LFP) L (IAN) _ y gon 
pnl [1] 


es decir, si existe una buena aproximación lineal de fen x, en el sentido ya conocido. 


Descompongamos el límite anterior según las funciones coordenadas de f, considerando por 


ejemplo f, , obtenemos 
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la función f, es diferenciable; y si la matriz asociada a A es A=Ía,), 18 [L2,---,m), 
jeX1,2,-"-,ny, la diferencial de f es 4,(h)=a,h, +a,,h,+..+a,,h,, lo que significa 
Y, 


d 
(recordando la expresión de la diferencial de una función real) que a,, = (x), 4), = e, a 


Ox, Ox, 


_of, 
, 4, dee 


n 


matriz Á no es otra que la matriz jacobiana de f, y la diferencial de fes 


(x). Y lo mismo para las restantes funciones coordenadas f,, f,,..., f,,. Luego la 


df (x)=Df(x)-h, 


o bien, llamando y = Fx), h=dx, 


dy, dx, 
dy, -Df a 
dy, dx, 


Es decir, la matriz de derivadas es localmente a las funciones diferenciables lo que la matriz 
(asociada a las bases canónicas) es globalmente a las aplicaciones lineales (una aplicación lineal es 
una función diferenciable cuya matriz jacobiana es constante). 

Composición de funciones diferenciables 


En este tema tenemos un teorema fundamental para todo el cálculo diferencial en varias variables. 


A A A $ 11 - * ea 2 
Teorema. Sean f y g dos funciones diferenciables (o de clase C*,o C”, re N ) definidas según el 
siguiente esquema: 


R n É R m g R P 


(x,, Le x,) ely, 3 a! plz, Sec 2,); 
entonces, la función compuesta go f : R"—> R” es diferenciable (o de clase C*, o C”). 
Aceptado este importante resultado podemos razonar como sigue. Sabemos que 


dz, dy, dy, dx, 


de donde 
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dz, dx, 
=DgDfA | 
¡de dx 


1) n 


pero, por la definición misma de diferencial de una función, sabemos que ha de ser, para la función 
diferenciable go f, 


dz, dx, 


luego 


que es la llamada regla de la cadena para las derivadas parciales de funciones compuestas de 
funciones diferenciables; resultando entonces además que 


d(go f)=dg0df: 
la diferencial de la función compuesta es la compuesta de las diferenciales. 


Es decir, en definitiva lo que tenemos aquí es una generalización de la propiedad (allí global) que ya 
conocíamos para las aplicaciones lineales (la matriz de la aplicación compuesta es el producto de 
las matrices), al caso (aquí local) de esas aplicaciones lineales particulares que son las diferenciales 
de las funciones diferenciables. 


Si recordamos la definición de las matrices jacobianas, 


a , ¡e ÁL2.... ph 

rl : E) Dígo f)= Es s je [1,2,...,n) 
Ox dx 

17 ke dl2..m) 


el 


el desarrollo de la matriz D[go A, de p filas y n columnas, nos da las derivadas parciales de la 


función compuesta, es decir, las de las funciones coordenadas Zz, respecto de cada una de las 
variables independientes Xx ;: 


dz, xx 0z, 0), 


o, 1 0Y, dx, 


Ejemplo 1. Supongamos que una función diferenciable, z = z(x, y), verifica en todo su dominio de 


definición la condición 


xz,-yz,=0, 
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y queremos saber en qué se traducirá esta condición haciendo el cambio de variables 


x=rcos0 


y=rsen0 
(coordenadas polares). Será 


] coso —rsenó 
X ) í 


L % 
tr adela 2? Pla, = 
FE “2 | send rcos0 
E (cosO z, +sen0z, —rsen0z,+rcosÓOz, ) ; 


y resolviendo el sistema, naturalmente /ineal (diferenciar es linealizar) en Z, y Z,, resulta 


z, sen0 
Zy rcosO rcosO z, —sen0 z, 
Z A O > 
] coso senó0 r 


=rsen09 rcosO 


coso 2, 
—=rsen0 2, rsen0z,+c0s0 z, 
ZE E A KKÁKÁKáKÁKÁ Y AA AOAOAOKÑ 
coso send r 


=rsen0  rcos0 

valores que llevados a la ecuación dada la convierten en 
Z¿ =0, 

que es la nueva forma que toma la condición. 
(Observación al margen. Esta nueva condición nos dice que z no depende de 0, sino sólo de Y, o 
dicho de otro modo, la gráfica de la función es una superficie de revolución de eje OZ; y puesto que 
el recíproco es inmediato, resulta que la condición dada es precisamente la característica diferencial 
de dichas superficies.) 


Ejemplo 2. Supongamos que tenemos la superficie 


Xx =UCOS V 
y=(1-u)jsenv  (u,v)e Rx[0,277), 


z=2u-1 
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: T 
y queremos calcular el valor de z, en el punto imagen de (uv) = 77 , es decir, en el punto 


(x, y, z) = 3, > 3 |. Tenemos: 


Dz Dz D(x,y) ] COSV —usenv 
A E A FA Z 
Y y 


-senv  (1-u)cosv 


a 
9 
SN 
5 
! 
O 
pe 
ÉS 
E 
== 
1! 
) 
a 
7 
ha 
_A 
a, 
ra 
BS 
= 
_— 
11 


o bien, 
Z, =COSVZ,—SEnvz, 


z, =-usenvz, +(l-u)cosvz, 


v 


sistema de Cramer del que obtenemos 


Z, —senv 
2 lz, (l-u)cosv|  (lI-ujcosvz, +senvz, 
] COS V sen v cos” v—u 


-usenv  (l-u)cosv 


y como z,=2, z,=0, el valor de la 


2(1=u)cos v 


derivada pedida es z = ; SN 
cos” vu 
en el punto en cuestión, 
1) 
EE 17 
HERA 2 AB 43 
Anas JA EN E 
TIRAN => h 
aa y E 


La figura muestra la superficie, cortada 


l 
por el plano y = 7 y la recta tangente 


a ella en el punto considerado y paralela a 
OXZ, es decir, contenida en dicho plano, 
recta que es también, naturalmente, 
tangente a la curva intersección y cuya 


443 


pendiente es z, =——=1,39, el valor 


5 


que hemos calculado. 


También podríamos haber utilizado 
directamente las diferenciales: 
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dx =cos vdu —usenv dv 
dy = —senvdu +(1—u)cos vdv 
dz =2 du 


De las dos primeras ecuaciones (sistema de Cramer) resulta 
dx  —usenv 


dy (1-u)cosv (1—u)cos v dx + usen v dy 


1 


du 


, 


Cos V —usenv cos” v—u 


-senv  (1-u)cosv 
lo que llevado a la tercera da 


, AL=u)cos v dx + usen y dy " 21—u)eosy y SE EN 


2 
cos” v—u COS VU cos” Y —u 


das 


obteniendo a la vez las dos derivadas parciales Z, y Z,. Particularizando en el punto fijado, 


, 


4/3 
= —— qx. 


dz 


' ; l 
De aquí tenemos inmediatamente, por ejemplo, que en ese punto (x, 2)=[43 E 3 |el 


gradiente de z (como función de x y de y) y el plano tangente son, respectivamente, 
44/3 443 
valtB ol 2 28-45), o bien 4/3x—52+3=0. 


Teorema de la función inversa 


Es análogo al que conocemos para funciones reales de una variable, haciendo aquí el jacobiano, o 
determinante de la matriz de derivadas, el papel que allí hacía la derivada. 


= ys v1 
Teorema. Sea f una función de R" en R”, de clase (al menos) C” y tal que en un punto 
a= la, ¡Borjas a,) el jacobiano, es decir el determinante de la matriz jacobiana, 


Le (a) = 7 (a) = det Df (a), 


no es nulo. Entonces existe en un entorno del punto h=(D,,b,,....b,)= f(a) la función inversa 


f", de clase (al menos) C*. 


Lo admitiremos sin demostración. 
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Como consecuencia de lo visto para la composición de funciones, resultará que, en un entorno de », 
"—| . =] . . ., . ze . . 
Df -Df" = D(/ of ) = Di, =1, donde ¡, es la función idéntica en R” e / su matriz 


jacobiana, es decir, la matriz unidad; luego 
Df” =(Df)": 


. . . ., . —| . . . . 
la matriz jacobiana de la función inversa f ” es la inversa de la matriz jacobiana de f. Y como 


consecuencia, también el jacobiano J pa es el inverso de /,. 


El teorema, que hemos visto en forma local, puede también hacerse global como en una variable. Si 
el jacobiano de f'no es nulo en ningún punto, fes biyectiva. En los ejemplos del punto anterior ya 
hemos topado con esto, aun sin decirlo. Cuando hemos dicho que los sistemas lineales que hemos 
encontrado eran de Cramer, estábamos implicando que las funciones que constituyen los cambios 
de variable son biyectivas. Concretamente, en el cambio a polares, 


p:R?*=>R? 
p(r,0)=(x, y)= (rcos0,rsen0), 


el jacobiano es, como hemos visto, /,, =r, y por lo tanto es biyectiva salvo en el origen; 


efectivamente, (0,0) = p(0,0) para todo O. 


En el otro ejemplo, el de la superficie en paramétricas, la aplicación (u, v) > (x, y,z) tiene como 


A ; 2 ; e Ñ 
jacobiano cos*v=u, luego es con seguridad biyectiva en todos los puntos en los que esta 
expresión no se anula, que en la superficie son todos los que no pertenecen a la curva 


_ o! e 3 , 1) 
(x,y,z) = (cos v, sen y, cos2v). 


Teorema de la función implícita 


Supongamos que tenemos una ecuación 
Pe (x, y) =0; 


está claro que esta ecuación establecerá una relación entre las variables x e y, es decir, podrá definir 
una de ellas como función de la otra. Gráficamente podríamos verlo como la intersección de la 


gráfica de la función z = F(x, y) con el plano z=0, lo que podría ser a su vez la gráfica, dentro 
de este último plano, de una o varias funciones y= p(x) Ó x= w(y) (o más simplemente 


y= y(x) Óx= x( y) ). De este problema se ocupa el llamado teorema de la función implicita. 


Teorema. Sea F una función de clase C” (m>1) en un entorno Y de un punto (a,b) e R?, con 
Fla,b)=0 y P, (a,b) + 0; entonces existe una función (, definida y de clase C” en un entorno 


Ude ae R, tal que p(a) =D y Ex, p(x)) = 0 para todo xe U . 
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Se dice entonces que esta función Y está definida 
implícitamente por la ecuación dada. 


Observación. Se habla también en este caso, 
impropiamente, de la función implícita (. 
Impropiamente, porque lo que caracteriza a esta 
función es sólo una manera de ser definida, no una 
naturaleza distinta de la de las demás funciones. 
Pero la costumbre ha impuesto la expresión. 


Ejemplo 1. Sea la ecuación 


[] 
IIA vs : A 
Ar F(x, y) = y —xy-3x+9=0, 
!] 


pipi 


que se verifica en el punto (x, y) 5 (3,3); 


F,(x,y)=3y*-x", F 


y 


(3,3) =18%0, 


es decir, se cumple la condición de existencia: la ecuación define en efecto, en un entorno de x=3, 
una función y = y(x) de clase C” como F y tal que y(3)= 3 (la figura representa la superficie 


z= F(x, y) , el plano z=0 y la curva intersección). 


Hay que subrayar desde ahora el carácter local de esta definición implícita 
de funciones. Así, en el caso del ejemplo, para x = 3 la ecuación 


F(3, y)= y* -9y=0 


admite otras dos raíces, y =0 e y=-3, es decir, la ecuación se verifica 
también en los puntos (3,0) y (3,-3), cumpliéndose en ambos la condición 


de existencia; luego la ecuación define implícitamente otras dos funciones 
en un entorno del mismo valor de x, como se muestra en esta segunda 
figura. 


Veamos otro ejemplo en el que la condición no se cumple. 


Ejemplo 2. Consideremos ahora la ecuación 


F(x, y)=6y'+x"y-17xy+18=0, 


que se verifica en (x, y) = (3,1) ; pero 


O: 


1 


F,(x,y)=24y"+x"-17x, F,(3,1) 
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la condición no se verifica. La figura muestra lo que 
ocurre en este caso. En cualquier entorno del punto 


(3,1) [rodeado de puntos] la curva intersección de la 


superficie y el plano no contiene ningún punto a la 
derecha y contiene dos puntos a la izquierda, por lo 
que no queda definida en él ninguna función 


y=p(x). 


También hay que observar que la condición del 
teorema no sólo se refiere a la existencia de la función 
implícitamente definida, sino también a su regularidad 


igual a la de F. Por ejemplo, la ecuación F (x, y)= x* +y*-1=0 se verifica en el punto 


(x, y)= (1,0), no cumple en él la condición del teorema puesto que F, (10)=0, y sin 


1 
embargo es evidente que define implícitamente la función y = li yt ) . Lo que ocurre es que, 


siendo F de clase C”, esta función no es ni siquiera derivable en el punto dado. 


Derivación de funciones implícitas 


Si F(x, y) = (0 define y como función de x, y = y(x), siendo la función compuesta E(x, y(x)) 


idénticamente nula, según la regla de la cadena podemos escribir 


x ) 


pp. Pe2)) - (A, E | Jona, 
y | 


o particularizando para el punto considerado, 
F, (a,b)+F, (a,b) y'(a)=0, 


lo que nos dará el valor de y (a) . Y lo mismo para derivadas sucesivas: derivando + Fy' =10 


obtenemos, 


E + Fs, y «+ $, y" => 0 » 


de donde obtendríamos el valor de y”, etc. 
Ejemplo. Volviendo a considerar la ecuación 


y =x*y-3x+9=0 


resulta, derivando y particularizando en el punto (8,3), 


3y y 2xy=x*y'3=0, 27y'-18-9y-3=0, y'(B)= 


Es 


derivando de nuevo de la misma forma y dando valores a x, ye y” resulta 


6yy?+3y?y"—2y-4xy'—x?*y"= 0, y"(3)= - ' 


etc. Con esto podríamos ya, por ejemplo, escribir y(x)=3+ = e 3) ss z (x— 3) + o(x e. 3) 


Caso de varias variables independientes 
Si la ecuación fuera 
FM Xy" % 399) Sl da 


tendríamos un teorema de existencia análogo; si la ecuación se verifica en un punto 
(a,,a,,-*-,4,,,b), con F, (a,,a,**,4,.,b)+%0, entonces queda definida una función 


O | 


PE YAA dr 1) en un entorno de (a,, a), a IA de la misma clase que F' y con 


a,,4,,"",0, ,)=b. Y en cuanto a la derivación, tendríamos análogamente que, por ser la 
y 1 2 n-1 5 q p 
función EA AA] de X,,X,,..., Xx, , idénticamente nula, 
1 0 0 
pr? Coto 9) _ (y 0. F,, E) e A 
a 1 dí n—1 y PS 
Dl ota) 0 0 1 
Ya, Yo, Ye, 
=(E+By, E+Fy, +“ FE +E)y, )=(0 0 0), 
lo que nos daría los valores de pd ; Ea di e 2 y nos permitiría seguir derivando. 


Ox, 


n—1 


Ox, 0x, 


(Observación. Por claridad de exposición hemos dado distinto nombre, y, a la variable cuya 
definición implícita como función de las demás queríamos estudiar. Pero, naturalmente, el 
problema podría plantearse igual eligiendo cualquier otra variable como presunta función, y no 
necesariamente la última.) 


Ejemplo. Sea la ecuación 


F(x,y,2)=2x* -3y?-2 +3y2+y+2=0, 
que se verifica en los puntos (11,1) y (11,2); como 


F.(x,y,2)=322+3y, F.(LL-1)=-3+3=0, F.(11,2)=-12+3=-9%0, 
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la ecuación define con seguridad, en un entorno de (x, y) = ALL), una función z = z(x, y), de 


clase C* como F, tal que z(1,1) = 2. Derivando respecto de x y particularizando en ( 1, 1,2) ' 


> 4 
4x-32"7,+3yz,=0, E, PE 
E: ao > 9 
y análogamente respecto de y, 
2 -0y +37 +1 l 
69-377, +32 4d, +1 =0, A, z,(1,1)=-. 
; 3z"-—3y : 9 


La figura muestra la gráfica de z y de su plano 
4x+ y +13 


ps 


tangente en el punto, z= 


Calculando directamente con diferenciales en lugar 
de derivadas parciales tendríamos 


Ñ 
A 4xdx — 6 ydy —-3z*dz +3 ydz + 3z2dy +dy=0, 
INNATO 
Is 
e 01/% 4x <6y +32 +1 
A A de = —_—_— MH AA BY, 
$ E E 32" =3y 
le y particularizando en (1,1,2), 
2 


dz di, 
g 9 


lo cual, recordando que dz=z,dx+z dy, nos da a la vez las dos derivadas parciales (y, 


sustituyendo dx, dy, dz por x=1l, y=1 y z-2 respectivamente, el plano tangente 


4x+ y+13 
z = —__—__—— 
9 


En cuanto a las derivadas parciales segundas, derivando otra vez respecto de x la primera expresión 
obtenemos 
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2 2 l 
4622. -327,,+3y2,,=0, 4- € a, =B, ZA) =—; 
: 3 81 % dl 243 


y lo mismo para z,, y Z,,. Esto nos permitiría, como más arriba, escribir el desarrollo de Taylor de 

la función z en dicho punto (lo que no es posible en este caso es hacer una representación gráfica 
: : 4 

del conjunto, puesto que nos estamos moviendo en KK” ). 
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Sistemas de funciones implícitas 


Supongamos ahora que la función /" que define la ecuación no es real, sino vectorial, es decir, de 
R" en R”,con n>M, pongamos n=M+ p: 


a ia eS m. 
Ple 0%,,, Jm )=0€ R > 
esto equivale, como sabemos, a tener un sistema formado por las m funciones coordenadas: 


A 
A 


El teorema de existencia es entonces análogo a los anteriores; si F está definida en un entorno Y de 
un punto (a,,...a,,b;,...b,, ) € R”, anulándose en él, y el jacobiano del sistema respecto de las 


presuntas funciones no es nulo en dicho punto, 


A 
) 


entonces existe definida por el sistema, en un entorno U de (a,--,a, Je R” una función de R? 


en R”, 
Je ylx,:=",x,) 
(es decir m funciones de R” en R, 
Yi =Yr (e,,-=,x,), k= 1,2,--,m), 


de la misma clase C” que F, con ylx,*-,x,)=b (o bien y, (x,,-,x,)= b,, k=1,2,-**,m), y 


tales que la ecuación dada, Elx, A OS 2 ) =0 € R” (0 bien el sistema 


Ed 2, (e) 1,)0 9, (0,0 x,))=0, J=12,--:,m) 


se verifica idénticamente en U. 


(Cabe hacer aquí la misma observación que más arriba en cuanto a la elección, y notación, de las 
variables presuntas funciones de las otras.) 


En cuanto al cálculo de las derivadas, se tratará otra vez de aplicar la regla de la cadena. 


Ejemplo. Sea el sistema 


2x* -8y-u*+v*=0 
3x0y+y*-u +2? =0 
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que se verifica en el punto (x, y,u, y) = (2,1,3,1) . Veamos si define dos funciones ul x, y) : v(x, y) 


en un entorno de dicho punto. 


El jacobiano es, en un punto genérico, 


y particularizando en (2, E 5 ; 


-6 53 
6 =-36+135=99%0; 


es decir, se verifica la condición de existencia. Derivando respecto de x, 


O 6x? —-2uu, +5v'v, =0 
9x”y-3u%u, +6vv,=0, 


sistema que, por ser su determinante el jacobiano, que no se anula en el punto considerado (ni, por 


continuidad, en un entorno del mismo), es de Cramer, y nos da 


bx sy 24 5 
-9x*y 6y? -36 6] -144+180 36 4 
Li E ¿AAA 
24  5v 6 S| -36+135 99 11 
—3u? 6y? —27 6 
Qu 6 —6 -—24 
R 34 -9x y (1.2 -36|_ 216-648 -432_ 48 
a Sn 99 99 99 11 
34? 6” 


Derivando ahora el sistema respecto de y, 


8-2uu, +5v*v, =0 
(1D) 


3 E 2 o 
dx" +4y"—3u'u, +6vv, =0, 
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8 Sy* 8 5 
_[3x-4y" 6v 5 (2 ¿Es ios 
d O 99 99 99 ” 
-3u* 6v? 
—2u 8 —6 8 
E 34% -—3x*-4y »,(2,1)= 18 al _168+144 _312_ 104 
d la mm + > 99 99 99 33 
34 6y? 


También habríamos podido operar diferenciando el sistema, proceso análogo pero más simple y 
sintético: 


— 2udu + Sv* dv =-6x*dx + 8dy 
— 3u* du + 6v*dv =-9x? ydx — (3x* +4y* )dy ¿ 


— 6x*dx + 8dy Sy? 6 $ 8 5y* 
91? ydx —(3x* +4 y*)dy 67 |-9x?%y 6v? -3x9-4y* 6y? 
du= dia. ENDE "O A OR d dy, 
—-2u  5v =2u  5Sv =2u  5v 
3 y —-3u? 6v? -3u* 6y? 


lo cual, puesto que du =u,dx +u,dy, nos da de una vez las dos derivadas parciales de u; y lo 


mismo si despejamos dv. 


En todo caso, derivando las expresiones de 4, y u,, o bien simplemente derivando de nuevo los 


sistemas (1) y (ID, obtendremos las derivadas segundas. Derivando (1) respecto de x, 


; 2x—2u ? —2uu,, +20vw 2 +5v*y, =0 


18xy—6u4 ? -3u%0,, +12vv?+6vv,, =0, 


y dando valores, 


o s= 0 [60 +50, = E 
121 + TZ] ú ió 121 
3 
6 E =0 Di, tg 


121 Sl S 121 
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e , le 48952 

31716 6 150764 1127 31716 

untaye  L0TS Ól_ 150764 , (pp 1]21 31716|_ 168000, 

$ 11 46 3 3993 ? 121 99 1331 
27 6 


y del mismo modo, derivando (1) y (ID) respecto de y obtendríamos las restantes cuatro derivadas 
segundas. Y así sucesivamente para derivadas de orden mayor. 
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lll — Fórmula de Taylor. Extremos 


La fórmula de Taylor proporciona, lo mismo que en funciones de una variable, una 
aproximación polinómica de una función, aproximación que, si la función tiene suficiente 
regularidad, será de grado mayor que la unidad, es decir, más allá de la aproximación dada 
por la diferencial. 


Fórmula de Taylor 


Sea f: D->RK una función de clase C” y con derivadas parciales hasta el orden n+1 en un 
dominio de /R?, D, que contiene en su interior el segmento ((x +2, yo +kt)[r e[0,1]), 


segmento que es la imagen de la función 
g: [0,1] >R 
ele) = (x(0).y(0) = (2, +», y, +10), 


de clase C” por estar definida por dos polinomios en la variable f. Sea ahora la función 


p=f>0g:[01]>R 
p(0)=(2(0) =$ (%, +ht, y, +M), 


que será derivable hasta el orden n+1 por serlo f. Apliquemos a y la fórmula de Taylor: 


OOOO o 


olt)= plo) + LL 
Pero 


o()= f(x, +h,y, +K), p(0)= F(x0,y0), 


y según la regla de la cadena, matricialmente 


Os 2. =(£,(x, +ht,y, +k0) £,(% +ht, 3,+40)| y )- 


=h f,(x, + ht, y, +) +k f, (2, + ht, y, +10), 


y para 1 =0, 


39 


p'(0) =»f, (xo, Yo) + M4, (%o, Yo) 
(que no es otra cosa que la diferencial de f). 


Derivemos ahora q” ; como 


. ls h 
Ll 4.1, +) Fo + ht, y, + kt) (> 


= hf o, (Xp + HL, yy +0) + Kf, (xp + ht, y, +Kt), 


y análogamente 
df, kt 
a O = Af (o + ht, Yo +4) + M5, (%o + ht, yy +1), 


resulta 


p"(1) =1f.(x, + ht, y, + kt) + 2hKf,, (x, + ht, y, + kt) +? f,, (xp + ht, yy + Kit), 


p"(0) =» f., (xo, Yo ) + 2AKf.,, (o, yo) + f, (xo, Yo). 


Esta expresión tiene una semejanza formal con el binomio de Newton: utilizaremos por ello la 
notación 


” () 
p"(0)=| 2f, (o, Yo) +4, (0.90) |” 
donde el superíndice indica potencias de » y k y órdenes de derivación de f. 


De forma análoga pero más farragosa iríamos obteniendo, con la misma notación, 


p”'(0) A A (X.Y) + ké, (0.3) =M Fo (0, Y.) +31Kf (Vo Yo) + 
+3hk" f,, (xo, Yo) HE e (o, Yo), 


y en general 


n (») 
y! (0) =[ Af, (o, Y) + Mf, (o, Yo ) | , 
expresiones a las que podríamos llamar potencias formales de la diferencial(*). 


Todo lo cual, llevado a la primera expresión, nos da el contenido - no enunciado hasta ahora - 
del teorema, es decir, la fórmula de Taylor para funciones reales de dos variables: 


En +h, y, + k) = Po») +, (20,90) +4, (20, Y ) | + 


1 : 1 s 
+5, (xo. Y.) + Ef, (%.-1)] ) ++ [1 (0,0) + 1, (2%: yo )] de 


(*) En realidad son las diferenciales de orden superior, d” f ; pero una introducción seria de esta noción nos llevaría 


demasiado lejos. 
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ll (20,90) +4, (20,10) ]7 ++2[1 (20,90) -+ 4, (30.90) > A 


bf Ax +0h y, +0K)+kf, (x,+0h, y, +0k)|">, 0<0<1. 


“a + ay ¡sl 
El último término, es decir, el llamado resto de orden n, r, (a, k), puede también escribirse 


1, (Mk) =| 


"s(h,k), con lim e(h, k)=0 


o bien, con notación de Landau, 


1. (1) =0 (Jo) 0 (17 +02. 


Naturalmente, si se tratara de funciones de un mayor número de variables obtendríamos una 
fórmula análoga; los términos, para tres variables por ejemplo, serían ahora 


1 m 
(0:70:20) +4, (o or20) +1. (or 3o:20) 7 con m=1,2,--,n 


Lo mismo que en funciones de una variable, el desarrollo de Taylor de grado n= 0, es decir, 
para dos variables, 


F Qro +h y, +k)- F(x0, 90) =hf, (x, +0h, y, +0k) + Kf, (x, +0h, y, +0k), 0<0<1, 


es la llamada fórmula de los incrementos finitos. Y lo mismo para más variables. 


Ejemplo. Sea la función 
2 


Xx 
,y)= arcí, , 
16) 


de clase C” en cualquier dominio del plano que no contenga puntos de la recta y =-2. 


Calculemos su desarrollo de Taylor de orden tres en el origen (0,0). 


2x(y+2) 
S(x,y)= po $, (0,0)=0 
(y +2) +(x?-1) 
241] 1 
AA E) 
: (y+2) (x? -1) 5 
y del mismo modo se va obteniendo 
1. (0.0)=<, f., (0,0) = 0, Fa k0,0)==>5 DE 5? Falo, == <> fonlo, ==, 


siendo nulas las otras dos derivadas de tercer orden; luego 
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Fx ) marche 3 eE E e 28 A + 
:, 2 Ms ans as alas +15 
y 24 2 hp 11 


+ olle A A A +o(x? +] 
: 2353 35 3 »>3 50375 j 


Extremos locales de funciones reales de varias variables 


Lo mismo que en el caso de funciones de una variable, una aplicación importante de la fórmula 
de Taylor es el estudio de los extremos locales de una función. Lo veremos a continuación para 
funciones de dos variables reales, generalizándolo después. 


Caso de dos variables 

Definición. Se dice que una función f definida en un dominio D =.R?, tiene un máximo 
[respectivamente mínimo] local en un punto (o, 7) € D si existe en D un entorno de dicho 
punto, V, tal que para todo (x, y)e V, con (x, y) 4 (xo, yo), se verifica la desigualdad 
E (x, y) <P (eo, yo) [respectivamente f (x, y) > f EM yo)1. 


Teorema. Para que una función que admite derivadas parciales en un punto tenga en él un 
extremo local, es condición necesaria que aquéllas se anulen en él (que el gradiente sea nulo). 
Es inmediato recordando la definición de las derivadas parciales como derivadas de las 


funciones de una variable f (x, y») y f (x, , y) , que tendrán un extremo si lo tiene f. 


Pero esta condición necesaria no es suficiente; basta pensar, por ejemplo, en la función 
E (x, y) = xy: sus derivadas parciales se anulan en (0,0), pero no tiene extremo en dicho 


punto; su valor en él es O y toma valores positivos en el primer y tercer cuadrantes y negativos 
en el segundo y cuarto. 


La fórmula de Taylor nos permitirá encontrar condiciones suficientes. Sea f de clase C ” en un 
dominio D, y (x,,y,)€ D un punto tal que f, (xo, Y,)=f,(%,, Y.) =0. El desarrollo de 


Taylor de orden uno en dicho punto será 
e 1,2, 
f(x, +h, yo +k)= (xo, y) +1 1.2, +0h,y, +0k)+ 
+2MH,, (x, +0h, y, +0k) +1? 4, (x, +0h, y, +0k) |. 


Y si la suma entre corchetes no se anula, su signo será, por la continuidad de las derivadas 
segundas, el mismo de 


VA. (xo, Yo) + 2hKf, (%o, Yo ) +, (%o, Yo), 


o matricialmente 


La (%o, Yo) Leg (0, Yo) 


forma cuadrática cuya clasificación haremos a la vista de su matriz. Llamando 


(» e PRI pela) 
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Lula) E A3d) (2 ha 
Paiñ FL At Ja) 4 4, 


el determinante de esta matriz, al que se llama el hessiano de fen (xo, y» ) , es 


H(%o,Y0)= 


2 
= 414 4), 


41 4) 
a 


A, 
resultando: 


1) si A (x,, Yo ) > 0 la forma cuadrática es definida, es decir, su signo no cambia al variar h y 
k, y es el mismo de a,, (y de a,,), luego hay un máximo local si a,, <Ú y un mínimo si 
En >U; 

2) si El (o, Yo ) < 0 la forma es indefinida, es decir, su signo cambia según los valores de k y 
k, y no hay extremo; 

3) si H lo. Yo ) = 0, la forma es semidefinida, es decir, se anula al menos para algunos valores 


de h y k, y como el valor de la forma cuadrática no es exactamente la variación de la 
función, sino que hemos despreciado un infinitésimo de orden superior al tomar las 


derivadas en (0, Y) y no en (x, +0h, y, + 0k) , es un caso dudoso: la variación de la 
función podría conservar o no su signo al variar h y k. 


Ejemplo. Sea la función 
F:R=>R 
fx, y) = seníx + y)+ cos(x - y). 


Busquemos los puntos donde sus derivadas parciales 
se anulan (puntos estacionarios): 


sn + y)- sen(x—y)=0 
cos(x + y)+ sen(x » y) =0 


es decir, sumando y restando, 


o x+y=2+ ka 


sen(x-—y)=0 yor 
r k-+k 
e A TT Tr 
4 2 y los puntos estacionarios son E K o —+ 12) ; 
Tk, —k, 4 24 2 
=—+ Hr 
4 2 


donde K y L son enteros de la misma paridad. El hessiano es, en un punto cualquiera, 
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H(x, y) y > sen(x Sud y)- cos(x = y) SS sen(x E y)+ cos(x Ek >) Ml 


—sen(x + y)+cos(x-y) — sen(x+ y)- cos(x-—y) 


= 4 sen(x + y)cos(x — y)= 2(sen2x + sen 2y), 


y en los puntos encontrados, 


az Af E + 12) — 2 se Kn) + sea + Lx) = e BEyE paa) 
4 24 2 2 . 4 (K y L impares), 


luego no hay extremo en el segundo caso, y en el primero tenemos, haciendo K =2k, L=21, 


—2 (% +1 par) 
p. (% +/ impar), 


| dr 7 
Í. ES ka, + in) 5 sen E +(4+ Dx) -cos((k -1)7)= | 


luego hay un máximo o un mínimo local respectivamente. 

Caso de tres variables 

Desde luego, tanto la definición como la condición necesaria de anulación de las derivadas 
parciales en el caso de que existan, son análogas cuando se trata de funciones de tres variables. 


En cuanto a la discusión, una vez localizados los puntos estacionarios EM Yo»2p) de la 
función, habrá que clasificar la forma cuadrática 


Ll dit) Loa) EAtBuktr) h 
(a k 1) Fo (%o» Yo» o) FLA Bor 2a) MEN k |, 
EPs) Lu Mio) EA l 


clasificación que dependerá de la matriz; llamando H a su determinante, el hessiano, y dando 
a los símbolos 4A,, y 4,, su sentido habitual, resulta: 


1) si AH>0, Aj, >0 y a,, >0, la forma cuadrática es definida positiva, y hay un mínimo 


local, 
2) si H<0, Ay >0 y a,, <0, la forma es definida negativa y hay un máximo local, 


3) en otro caso, con H 4 0, la forma es indefinida y no hay extremo; 

4) si H=0, la forma es singular, semidefinida, y es un caso dudoso (que se resuelve si las 
raíces no nulas de la ecuación característica tienen signos distintos, en cuyo caso la forma es 
indefinida y no hay extremo). 


Ejemplo. Sea la función 


F:R>R 


Fx, y,z) =x + y? +4x2+z* +4x. 
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3x? +47+4=0 ó á 6x 0 4 

2y=0 : puntos (2,0,-4) y (2,0,-2); H(x, y,z) = Ú Z QU: 

4x+2z=0 4 02 
1 04 

H(2,0,-4)=|0 2 0/=16>0, 4, = a y =24>0, a =12>0: 


4 0 2 


mínimo local en la, 0, — 4); 


0 4 40 
2 0|=-16<0, ta=Í: q =8>0, a41=4>0: 
o 2 


O 
no hay extremo en 3" $ 3 d 


Generalización a n variables. Los casos vistos de dos y tres variables no son más que casos particulares 
del general de n variables, que puede analizarse recordando algunos conceptos de álgebra. El término 
cuadrático del desarrollo de Taylor de una función de » variables es la forma cuadrática 


de matriz ¿- (a ) donde a. es el valor de la derivada f_. en el punto estacionario estudiado, matriz 
y dd 
que es diagonalizable, es decir, mediante un cambio de ejes la forma se convierte en reducida, $ A ne a 
k=1 


donde los coeficientes A, son las raíces de la ecuación característica, det(4 - AL) =0; la forma es 


entonces indefinida si entre dichas raíces las hay positivas y negativas (no hay extremo), definida positiva 
si son todas positivas (mínimo), definida negativa si todas negativas (máximo) y semidefinida si alguna es 
nula y las restantes tienen el mismo signo (caso dudoso). Y esto se traduce en la matriz así: si 


a tm e *s 
E (Y) Ma rr, det A=H (hessiano) 


as e 


son todos estrictamente positivos, hay un mínimo; si son alternativamente negativos y positivos a partir de 
a,, <0, un máximo; en otro caso, con H 0, no hay extremo; y si H =0 habrá que investigar las raíces 


no nulas de la ecuación característica; si las hay de distinto signo, no hay extremo. Es la regla que hemos 
visto más arriba particularizada para los casos n=2 y n=3. 


Extremos condicionados: multiplicadores de Lagrange 


Supongamos que necesitamos conocer los extremos locales de una función real de varias 
variables, pero no en todo su dominio de definición, sino restringida a una parte M de éste 
definida por una o varias condiciones, o ecuaciones de enlace. Es éste el problema de los 
llamados extremos condicionados, para resolver el cual disponemos de un teorema, debido a 
Lagrange, que enunciamos primero en el caso particular más sencillo. 
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Teorema. Sea f una función real de clase C * en un dominio Dc R?, una función de D en 
R también de clase C', M= ((y)|o(z y) = 0] cD y (a,b) un punto de M en el que la 


restricción de fa M tiene un extremo local. Entonces existe un único número real 4 tal que, si 
es 


F=f+10, 


se verifica 


OF OF 
——(a,b)=——(a,b)=0. 
00) (a0)=0. 


Admitiremos este teorema sin demostración. 


En él se basa el método llamado de los multiplicadores de Lagrange, consistente en resolver el 
sistema formado por las tres ecuaciones 


F.(x,y,4)=0 
Br (1,y,4) =0 
p(x,y)=0 


con las tres incógnitas x, y, 4. La resolución de este sistema nos da los puntos en los que 
pueden existir extremos. Para saber si esto es así y de qué extremos se trata, es frecuente que 
baste con considerar las condiciones del problema. 

Ejemplo 1. Supongamos que sobre la curva 


xx +8y? -3xy-4=0 (-4<x53) 


se ha construido un muro cuya altura en 
cada punto es 


ho y)=44 045, 


y queremos saber cuál es su altura 
máxima. La función de Lagrange será 


Flo y, 2) 24404 y + ale +8y -3xy-4), 


de donde 


> +Abx? -3y)=0 ' ( 


de stos Mk as 
1+ 4(24y? -3x)=0 2Bx?-3y) 24y? -3x E sc 


x* +8y" -3xy-4=0 
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Xx 
2 4x* -3x?-8=0; 
x +8" -3xy-4=0 


(x-2y)(Qx+4y+1)=0 a y= 
xx +8y -3xy-4=0 


raíz real única, que numéricamente aproximada es 


x, =1,566, de donde y, =0,783 y Mx,,y,)=5,566 ; 
1 Y py 


_ 2x+l 
b) 4 no hay soluciones reales; 


x* +8y"* -3xy-4=0 


luego el máximo pedido es 5,566, valor que se alcanza en 


el punto E , y,) (nótese que lo que hemos resuelto es un 
problema de extremos globales). 


Enunciemos ya el teorema de Lagrange en su forma más general. 


Teorema. Sea f una función real de clase C* en un dominio D <= R”, ( una función de D en 
R" (m<n,) también de clase C* (siendo Oi, P>)»»»>Q,, sus funciones coordenadas), 
M= ¡(3108 Je (2 %p:5,2,) =0€ R”) =D y a=(a,a,,...,a,) un punto de Men 


el que la restricción de f a M tiene un extremo local (siendo la matriz jacobiana (mxn) 


D(0,,0,,...,P,,) de e ca 
== de rango máximo m para que las m condiciones 0; [ptas JE, 
Dio) 


j =1,2,...,m sean independientes). Entonces existen m números reales únicos A, A,,..., A, 


tales que, si es 
F= +40 +4,P), + +A Dm 


se verifica 


OF OF OF 

2 laa E (a) (a)=0. 
Xx Ox, Ox, 

En este caso, naturalmente, al aplicar el método 

tendremos un sistema de m-+m ecuaciones (las n 

derivadas y las m condiciones) con m=>+m 


incógnitas (las 1 coordenadas y las m A,). 


Ejemplo 2. Supongamos que se trata de calcular 
las distancias máxima y mínima de la curva 


(elipse) 


e +yoz-1=0 
x+y+22-2=0 


al origen de coordenadas. 
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La figura representa la elipse como intersección del hiperboloide de revolución de una hoja y 
del plano cuyas ecuaciones constituyen el sistema dado. 


Tomaremos como función no la distancia sino su cuadrado, función polinómica más sencilla y 
diferenciable en todos los puntos, y que alcanza sus extremos en los mismos puntos que 
aquélla. Es decir, la función cuyos extremos buscamos es 


Hay z=x+y +2 
restringida a los puntos de la curva, y la función de Lagrange será 


Flx,y,z A p)=x+y +2 +2(x? 4 yz -1)+ u(x+y+22-2); 


NN CA A 
2y+42y+u4=0 

2x-2+A(2x+2)=0 2x—2+A(2x+2)=0 
22-422+2u4=0 
AE +y-2-1=0 +yz-1=0 
x+y"-z -1=0 

x+y+22-2=0 x+y+22-2=0 
x+y+22-2=0 


RARE POR 
-2144+8 
a) 32x?-2?-1=0 31?+2x-2=0 x= ÁS 143; 


x+z-1=0 z=1l-x 
puntos estacionarios: (14 43, Tf, 23-343), (- 1- 13, | =413, 2+/3); 


A=-1 
z=0 z=0 
2x-z-(2x+z)=0 ni ” 
E x"+y"-1=0 32x —4x+3=0, 
x+y-z"-1=0 
x+y-2=0 y=2-x 
x+y+22-2=0 


que no tiene soluciones reales; los únicos puntos son los dos hallados, cuyas respectivas 
distancias al origen son 


d, =/24-24/3)+7-44/3 = (15-843 =10694, d, =/15+843 =5,3718: 


la primera es, pues, la mínima y la segunda la máxima. 
Aplicación al cálculo de extremos globales 


Sea f una función de dos variables definida en una región A del plano limitada por una curva 
pl, y) =0, con f y y de clase C' y A cerrada (es decir, que incluye su frontera OA formada 


por los puntos de la curva); existen entonces ambos extremos globales de fen A. Y es obvio que 
si ftoma su valor máximo (o mínimo) en un punto interior de 4, en ese punto tendrá también un 
extremo local y sus derivadas parciales se anularán, y si lo toma en un punto de la curva 
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frontera, tendrá en ese punto un extremo condicionado, con olx, y)= O como condición de 


enlace. Luego bastará con buscar los puntos que verifican esas condiciones (si hubiera algún 
punto aislado en el que f no tuviera derivadas o tuviera una sola, nula, habría que incluirlo 
también) y comparar los valores que la función toma en ellos. (Si la función — acotada — está 


extendida a todo RR? sólo habrá que considerar puntos de derivadas nulas o inexistentes; por 
ejemplo, los extremos de la función de la página 42 son max f (x,y)=2, min f (x, y) =-2). 


extendida a todo R? sólo habrá que considerar puntos de derivadas nulas o 
inexistentes; por ejemplo, los extremos de la función de la página 42 son 
max f(x,y)=2, min f(x, y)=-2). 
Ejemplo. Supongamos que el recinto 


A=((,y)|" + y -2x-2y <510) 


está cubierto por una bóveda cuya altura en 
cada punto es 


(longitudes en metros), y cerrado por un muro perimetral vertical, y necesitamos 
conocer la altura máxima de la bóveda y la máxima y mínima del muro. 


2x+12=0 
2y=0 
36-72 103 


E = nh Y para el muro (multiplicadores de Lagrange), 


F(x,y,4)=20-(x +)? +12x)+ 2(x* + y? -2x-2y-510); 


Tendremos, en el interior de 4, | es decir, el punto (-6,0), en el que la 


altura es z(-6,0)=20- 


HE 02) =0 
0 0 y x+6 
de am. 525. yá 
3 +4(y-2)=0 x—1 y-1 7 


Ped si=0 x*+y?-2x-2y-510=0 [1 +y*-2x-2y-510=0 


2 
pe PE y AO 100, 25x? -50x-12519=0, 
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50+/2500+1251900 5+112 E o) 107 >) 
Xx = ——————_ >= ; dos puntos, | —,—| y | -—,-— |, en los 
50 5 S 3 5 5 
que las alturas del muro son dE) E En resumen: 
5 5 10 5 5 30 
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altura máxima de la bóveda, o =20,60 m.; altura máxima del muro, pd =16,57 m.; 


y altura mínima del muro, = =35,90 m. 


Este método se puede aplicar también, teóricamente, a funciones de más variables. Para tres, 
por ejemplo, si tenemos una función f : 4A>R, siendo A el recinto cerrado limitado en R? 


por la superficie p(x, y,z) =(), y queremos calcular sus valores máximo y mínimo, bastará 
hallar los puntos estacionarios en A y los de 04 donde puede haber extremos (problema de 


Lagrange) y comparar los correspondientes valores de f. Es evidente que, en la práctica, el 
aumento del número de variables hará los cálculos progresivamente abrumadores. 


Ejemplo. En la página 43 estudiábamos la función f (x, y, z) =x* + y? 4412 +7? +4x, y 
llegábamos a la conclusión de que tenía un solo extremo local, un mínimo en (2, O, 4); 


calculemos sus extremos globales en la esfera E de centro el origen y radio 5, 
+ y + z? <25, en cuyo interior se encuentra el punto citado. Veamos la frontera: 


F(x,y,2,4)=x + y +4x24 22 +4x4 4 (x* + y? +2? 25); 
3x* +42+4+24x=0 
2y(1+4)=0, 4=-1 o' y=0 
4x+2z+21z=0 
*+y+2-25=0 


3? +42+4-2x=0 
a) A=-1: ¿x=0 


eo...» 


47+4=0 
y +2 =25 


3x? +42 +4+24x=0 
puntos (o, 246, -1) y (0, 2/6, -1); b) y=0: ¿4x+2z2+212=0 


x+z2?=25 


3x2 +42? +42-4x? -2xz=0 


_ 3x7 +42+4_4x+2z 
x z | 


x+20=25 


sistema que, resuelto numéricamente, da cuatro puntos más; con tres decimales, 


(4,818, 0, 1,339), (2,173, 0, -4,503), (-1,717, 0, -4,696) y (-4,885, 0, 1,065); 


valores de £ f (2,0,-4)=0, f(0,+246,-1)=25, f(4,818,0,1,339) =158,661, 
F (2,173,0,-4,503)=0,090, f (-1,717,0,-4,696) =42,374, f (-4,885,0,1,065) = 


—155,809; luego el máximo y el mínimo globales son, con tres decimales, los números 
158,661 y -155,809; o bien, por ser f continua, f (E) =[-155,809, 158,661]. 
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